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Abstract 

We consider linear and semi-linear stochastic partial differential equations and establish 
f*^ ' stability in rough path sense. The main example we have in mind is the Zakai equation from 

£vj , the theory of non-linear filtering; in this case our results can be viewed as "robustification", 

applicable even when the noise in the observation process depends of the signal. 

oS'. 

Ph ■ 1 Introduction 

+H ■ 1.1 Motivation from filtering 

Stochastic filtering is concerned with the estimation of the conditional law of a Markov process, 
given observations of some function of it. The classical formulation looks at the case where the 
process is of diffusion type and splits into a first part (known as the signal) and a second part, 
known as the observation process with values in a vector space, and whose martingale part has 
stationary increments independent of the signal. (This amounts to a = in what follows and we 
will speak of the uncorrelated case.) However, very often filtering problems present themselves with 
the noise in the observation process being dependent of the signal. In this case it is much harder to 
achieve robust and stable filtering as the transformation; we will speak of the correlated case. To 
be specific, let B, B be independent multidimensional Brownian motions and consider 

o ; 

o . d 

dX = V (X)dt + ^2v i (X)dB l + J2 (J j( x ) dBJ (signal) 

i j=l 

dZ = r/ (X) dt + dB e M d (observation). 

Under generic conditions, the (unnormalized) conditional density ut of X t given a (Z s : s < t) 
is then known to satisfy Zakai's equation 

d 

du t = G*u t dt + Y^ N*u t dZ{ 



[g* -\J2N*N;\u t dt + j^N;u t ozi 



1 The present discussion is purely motivational and we shall not give precise conditions on the coefficients. 



where G = (V ,D-) + (1/2) Tr [(VV T + aa T ) D 2 -] is the generator of X and Nju = (<tj (x) , Du) + 
rjj (x) u. In other words, u = u(t, x; uj) satisfies a stochastic partial differential equation (SPDE) of 
parabolic type, namely 

d 
du = L (x, u, Du, D 2 u) dt + V^ Aj (x, u, Du) o dZ^ . 

i=i 

Note that we have written the observation differentials in Stratonovich form; this will be convenient 
for the rough path approach taken later in this paper. 

Classical robustification: In the uncorrelated case (a = 0) one has Mj = Mj (x, u) — rjj (x) u 
and it is part of the folklore of the subject that Zakai's equation admits a pathwise solution. Indeed, 
take a smooth path z and solve the auxilary differential equation <fi = <p ^ ■ 77 (x) dz J = <f>r] • dz. 
The solution is given by 

4>t = ^0 ex P ( / v ( x ) - dz J =<t>o ex P (v 0) • z t) 

and induces the flow map <j)(t,4> ) := <f> exp (rj (x) ■ z t ); observe that these expressions can be 
extended by continuity to any continuous path z such as a typical realization of Z. (w). The 
point is that this transform allows to transform the Zakai SPDE into a random PDE, sometimes 
called the Zakai equation in its robust form: it suffices to introduce v via the "outer transform" 
u (t, X) = <f> (t, v (t, x)) which leads immediately to 

v (t, x) := exp (—ri (x) ■ z t ) u (t, x) . 

An elementary computation then shows that v solves a linear PDE, in the sense that ^L is affine 
linear in v, Dv, D 2 v with coefficients that will depend on z resp. Z t (w), 

dv = L [x, v, Dv, D v) dt. 

Moreover, one can conclude from this representation that u — u(z) is continuous with respect to the 
uniform metric |^ — -^loo-ro Tl* This provides a fully pathwise "robust" approach to Zakai's equation. 
Rough path robustification: The classical robustification does not work in presence of 
correlation (a ^ 0). In fact, we cannot expect PDE solutions to 

d 
du = L (x, u, Du, D u) dt + y Aj (x, u, Du) dz, 

3=1 

surely well-defined for smooth z : [0,T] -> R d , to depend continuously on z in uniform topology. 
Our main result is that u = u (z) is continuous with respect to rough path metrico That is, if 
(z n ) C C 1 ([0,T] ,R d ) is Cauchy in rough path metric then (u n ) will converge to a limit which ill 
be seen to depend only on the (rough path) limit of z n (and not on the approximating sequence). 
As a consequence, it is meaningful to replace z above by an abstract (geometric) rough path z and 
the analogue of Lyons' universal limit theorem |11) holds. Let us briefly mention that this leads to 
numerous (essentially immediate) corollaries for our class of SPDEs, including Stroock-Varadhan 



2 Two (smooth) paths z, z are close in rough path metric iff z is close to z AND sufBcienly many iterated integrals 
of z are close to those of z. More details are given later in this article as needed. 



type support theorems, Freidlin-Wentzell type large deviations, maximum principle and splitting 
results. Since this has been rather well understood in the (rough path) literature, we shall not go 
into explicit detail here. 

At last, we should mention that the relevance of rough path theory for filtering in presence of 
correlation appears as motivating conjecture in Lyons' seminal 1998 paper [TO] on rough paths, 
a theory that proposes to understand paths together with their iterated integrals ("areas") in 
order to guarantee robustness of differential equations driven by such paths. This conjecture is 
supported in particular by the works of I. Gyongy [5J [7] who obtained explicit results of how 
"correct" approximations to the observation process, but wrong approximation of its area process, 
lead to convergence to the "wrong" solution. A quantitative form of this conjecture is the assertion 
that solutions to the filtering problem (i.e. the afore-mentioned SPDEs) depend continuously on 
the observation process seen as a rough path, the proof of which is the content of this paper. 

1.2 Structure and outline 

Consider the SPDE 

d 

du + L (t, x, u, Du, D 2 u) dt = ^ A fc (t, x, u, Du) o dZ k , (1) 

fc=i 

with scalar initial data u (0, •) = uq (•), assumed to be BUC (bounded and uniformly continuous) on 
R™, and give meaning and rough path stability for the solution to the above equations as a function 
of Z, the (canonical semi-martingale) enhancement of the d-dimensional observation process Z . 
Here L is a linear second order operator of the form 

L (t, x, r,p, X) = -Tr [A (t, x) ■ X]+b(t,x) ■ p + c (t, x, r) . 

Actually, our approach allows to deal with L semi-linear in the sense that / above may depend 
non-linearly on r with no extra effort. The first order operators A k = Afc (t, x, r,p) are affine linear 
in r,p. That is, 

A fe (t, x, r,p) = (p ■ <7fc (x)) + r v k (t, x) + g k (t, x) . 

We shall prefer to write the rhs of ([JJ in the equivalent form 



d\ d 2 d 3 

i=i j=i /t=i 



(Du ■ a l (x)) o dY t hl + u^ Vj (t, x) o dY 2 '^ +^g k (*, x) o dY 3 < k 



where Y = (Y 1 , Y 2 , Y 3 ) is a (d\ + d^ + c^-dimensional semimartingale. Our approach is based on 
a pointwise (viscosity) interpretation of (JTJ: we successively transform away the noise terms such 
as to transform the SPDE, ultimately, into a random PDE. The big scheme of the paper is 

Transformation lli In ,i/ t^x • i- i irl i 

u i— y u where u has the (gradient) noise driven by Y removed; 

i Transformation 2i9i 12i ,1 -• • i- , t^2 i 

u h- > u where u has the remaining noise driven by Y removed; 

19 Transformation 3--, ~, .-, .. . n . , ^^ 1 

u t-> u where u has the remaining noise driven by Y removed. 

None of these transformations is new on its own. The first is an example of Kunita's stochastic 
characteristics method; the second is known as robustification (also know as Doss-Sussman trans- 
form); the third amounts to change u 12 additively by a random amount and has been used in 



virtually every SPDE context with additive noise|f| What is new is that the combined transforma- 
tion can be managed and is compatible with rough path convergence; for this we have to remove 
all probability from the problem: In fact, we will transform an RPDE (rough PDE) solution u 
into a classical PDE solution u in which the coefficients depend on various rough flows (i.e. the 
solution flows to rough differential equations) and their derivatives. Stability results of rough path 
theory and viscosity theory, first introduced for a class non-linear problems [2], then play together 
to yield the desired result. Upon using the canonical rough path lift of the observation process 
in this RPDE, a truly (rough)pathwise point of view, one has constructed a robust version of the 
SPDE solution in equation. We note that the viscosity/Stratonovich approach allows us to avoid 
any ellipticity assumption on L; we can even handle the fully degenerate first order case. In turn, 
we only obtain BUC solutions. Stronger assumptions would allow to discuss all this in a classical 
context (i.e. u would be a C 1,2 solution); SPDE solution can then be seen to have certain spatial 
regularity. 

In the presence of gradient (rough) noise only this approach was implemented in [1] . We should 
remark that the usual way to deal with ([I]), which goes back to Krylov, Rozovskii, Pardoux and 
others, is to find solutions in a suitable functional analytic setting; e.g. such that solutions evolve 
in suitable Sobolev spaces. Interestingly, there has been no success until now (despite the advances 
by Gubinelli-Tindel [6] and Teichmann |13j ) to include (J]) in a setting of abstract rough evolution 
equations on infinite-dimensional spaces. 

Acknowledgement 1 The first author would like to thank the organizers and participants of the 
Filtering Workshop in June 2010, part of the Newton Institute's SPDE program, where this work 
was first presented. The second author would like to express his gratitude for a Newton Institute 
Junior membership grant. 

2 Recalls on viscosity theory and rough paths 

Let us recall some basic ideas of (second order) viscosity theory [3j H] and rough path theory 
[HJ [12]. As for viscosity theory, consider a real- valued function u = u (t, x) with t € [0,T] , x G K n 
and assume u 6 C 2 is a classical subsolution, 

d t u + F (t, x, u, Du, D 2 u) < 0, 

where F is a (continuous) function, degenerate elliptic in the sense that F (t,x,r,p, A + B) < 
F (t, x, r, p, A) whenever B > in the sense of symmetric matrices. The idea is to consider a 
(smooth) test function p and look at a local maxima (t, x) of u — p. Basic calculus implies that 
Du (t, x) = Dip (i, x) , D 2 u (t, x) < Dip (t, x) and, from degenerate ellipticity, 

d t ip + F(t,x,u,Dip,D 2 p) <0. (2) 

This suggests to define a viscosity supersolution (at the point (x, i)) to dt + F = as a continuous 
function u with the property that ([2|) holds for any test function. Similarly, viscosity subsolutions 
are defined by reversing inequality in ((5); viscosity solutions are both super- and subsolutions. A 



•^Transformation 2 and 3 could actually be performed in 1 step; however, the separation leeds to a simpler analytic 
tractability of the transformed equations. 



different point of view is to note that u (t, x) < u (i, x) — <p (t, xj+cp (t, x) for (t, x) near (i, x) . A 
simple Taylor expansion then implies 

u (t, x) < u (t, x) + a (t — t) + p ■ (x — x) + — (x — x) ■ X ■ (x — x) + o f \x — x\ + \t — t\ ) (3) 

as \x — x\ + \i — t — > with a = dt<p (i,x) , p — Dip (t, x) , X = D 2 ip (t, x) . Moreover, if ^ holds 
for some (a,p, X) and u is differentiable, then a = dfU (t, X) , p — Du (t, x) , X < D 2 u (t, x), hence 
by degenerate ellipticity 

d t ip + F(i,x,u,p,X) < 0. 

Pushing this idea further leads to a definition of viscosity solutions based on a generalized notion 
of "(dtu, Du, D 2 u)" for nondiffcrcntiable u, the so-called parabolic semijets, and it is a simple 
exercise to show that both definitions are equivalent. The resulting theory (existence, uniqueness, 
stability, ...) is without doubt one of the most important recent developments in the field of partial 
differential equations. As a typical resullo, the initial value problem (<9 t + F) u = 0, u (0, ■) = 
u e BUC(M") has a unique solution in BUC([0,T] x R") provided F = F{t,x,u,Du,D 2 u) is 
continuous, degenerate elliptic, proper (i.e. increasing in the u variable) and satisfies a (well-known) 
technical conditional. In fact, uniqueness follows from a stronger property known as comparison: 
assume u (resp. v) is a supersolution (resp. subsolution) and uq > vq; then «>!ion [0,T] x R n . 
A key feature of viscosity theory is what workers in the field simply call stability properties. For 
instance, it is relatively straight-forward to study (dt + F) u = via a sequence of approximate 
problems, say (dt + F n ) u n — 0, provided F n — > F locally uniformly and some apriori information 
on the u n (e.g. locally uniform convergence, or locally uniform boundednesaj. Note the stark 
contrast to the classical theory where one has to control the actual derivatives of u n . 

The idea of stability is also central to rough path theory. Given a collection (V\, . . . ,Vd) of 
(sufficiently nice) vector fields on R" and z g C 1 ([0, T] ,R d ) one considers the (unique) solution y 
to the ordinary differential equation 



y 



[t) = Y,V l {y)z l {t) 1 y(0)=y £R n . (4) 



The question is, if the output signal y depends in a stable way on the driving signal z. The answer, 
of course, depends strongly on how to measure distance between input signals. If one uses the oo 
norm, so that the distance between driving signals z, z is given by \z — ^I^.tq t i, then the solution 
will in general not depend continuously on the input. 

Example 2 Take n = l,d = 2, V = (Vi, V2) = (sin (•) , cos (•)) and yo = 0. Obviously, 

z n (t)= [ - cos (2mi 2 t) , - sin (2nn 2 t) ] 
\n n J 

converges to in co-norm whereas the solutions to y n = V (y n ) z n ,yQ = 0, do not converge to zero 
(the solution to the limiting equation y — 0). 



4 BUC (...) denotes the space of bounded, uniformly continuous functions. 

5 (3.14) of the User's Guide [3]. 

6 What we have in mind here is the Barles-Perthame method of semi-relaxed limits [4]. 



If \z — ^ I oo-rg t] ^ S replaced by the (much) stronger distance 

\ Z ~ Z ll-var;[0,T] = SU P 2^ \ z U,t i+ i ~ Z ti ,t l+ l > 

it is elementary to see that now the solution map is continuous (in fact, locally Lipschitz); however, 
this continuity does not lend itself to push the meaning of (|4]): the closure of C 1 (or smooth) paths 
in variation is precisely W ' , the set of absolutely continuous paths (and thus still far from a 
typical Brownian path). Lyons' theory of rough paths exhibits an entire cascade of {p- variation or 
l/p-H61der type rough path) metrics, for each p e [l,oo), on path-space under which such ODE 
solutions are continuous (and even locally Lipschitz) functions of their driving signal. For instance, 
the "rough path" p- variation distance between two smooth K d - valued paths z, z is given by 



i/V 



/ v^ U) ~(i) 

max sup > U t V,, - z t V, , 

j=l,...,[p] \(U)C[0,T] I 



where z s t = J dz ri (&- ■ -®dz r . with integration over the j-dimensional simplex {s < n < • • • < rj < i). 
This allows to extend the very meaning of ((U), in a unique and continuous fashion, to driving sig- 
nals which live in the abstract completion of smooth Revalued paths (with respect to rough path 
p- variation or a similarly defined l/p-H61der metric). The space of so-called p-rough pathqj is 
precisely this abstract completion. In fact, this space can be realized as genuine path space, 



C" 



0,p-var A Qj T ] ) G [p] ^dj\ Tesp C 0,l/p-H61 fa y] > Q ]p] ( R tf 



where G^ (K d ) is the free step-[p] nilpotent group over R d , equipped with Carnot-Caratheodory 
metric; realized as a subset of 1 + t' p ' (R d ) where 

ti p l (R d ) = R d © (R d )® 2 © • • • © (R d f lp] 

is the natural statespace for (up to [p]) iterated integrals of a smooth Revalued path. For instance, 
almost every realization of d-dimensional Brownian motion B enhanced with its iterated stochastic 
integrals in the sense of Stratonovich, i.e. the matrix-valued process given by 

5(2) : = ( / B l odBi\ (5) 

V-A) J i,je{i,...,d} 

yields a path B (u>) in G 2 (R d ) with finite l/p-H61der (and hence finite p- variation) regularity, for 
any p > 2. (B is known as Brownian rough path.) We remark that B^> = ^B <g> B + A where 
A := Anti (B^ 2 *) is known as Levy's stochastic area; in other words B (to) is determined by (B, A), 
i.e. Brownian motion enhanced with Levy's area. A similar construction work when B is replaced 
by a generic multi-dimensional continuous semimartingalcs; sec [5, Chapter 14] and the references 
therein. 



7 In the strict terminology of rough path theory: geometric p-rough paths. 



3 Transformations 

3.1 Inner and outer transforms 

Throughout, F = F (t, X, r,p, X) is a continuous scalar-valued function on [0, T] xR™ xlxl" xS (n), 
§ (n) denotes the space of symmetric n x n-matrices, and F is assumed non-increasing in X 
(degenerate elliptic). Time derivatives of functions are denoted by upper dots, spatial deriva- 
tives (with respect to x) by D.D 2 , etc. Further, we use (., .) to denote tensor contraction^!, i.e. 

Lemma 3 (Inner Transform) Let z e C 1 ([0, T] , W l ) , a = (ai, . . . , <r d ) C C h 2 (R n , R n ) and ijj = 
ip{t,x) the ODE flow of dy — V (y) dz, i.e. 

d 

^(t,x) =^2<Ti(i}{t,x))zi, i>(0,x) =i£l". 
»=i 

Then u is a viscosity subsolution (always assumed BUC^ of 

d 

d t u + F (t, x, r, Du, D 2 u) ~^{Du- a, [x]) z\ = 0; u (0, .) = u Q (.) (6) 

iff w (t, x) :— u (t, tp (t, x)) is a viscosity subsolution of 

dtw + F^ (t,x,w,Dw,D 2 w) = 0; w (0, .) = u (.) (7) 

where 

F^ (t,x,r,p,X) = F(t,tp t (x),r, (p^D^ 1 ^^)) ,(X,Dij)t\^ x ) <8> D^l^x)) + (p.D 2 ^ 1 ^^))) 
and 

/ * h\ k—l,...,n . H u. k=l,...,n 

D^X- d{f \ { '' X)) ) -BV-" - ^ M)) ^ 



dx l j " ' rt lx \ dx'xi / 

V / i— l,...,n \ / i,j=l,...,n 

XTie sa?7ie statement holds if one replaces the word subsolution by supersolution throughout. 

Corollary 4 (Transformation 1) Let ip = ip(t,x) be the ODE flow of dy = V(y)dz, as above. 
Define L = L (t, x, r,p, X) by 

L = - Tr[A (t, x) ■ X]+b(t,x) -p + c (t, x, r) ; 

define also the transform 

L^ = -Tr [A* (t, x)-X}+ b^ (t, x) ■ p + c^ (t, x, r) 



'We also use • to denote contraction over only index or to denote matrix multiplication. 



where 

A*(t,z) = (A(t,^ t (a;)),-DV'r 1 U t (x)«'-DV't" 1 |^(x)> I 
h 4 ' (t,x) -p = b{t,ip t (x)) ■ {p,Dil>t\ t ( x) ) - Tr(A(t,ip t ) ■ (p, D 2 ipt\ t ( x ))) , 
c^{t,x,r) = c{t,ip t (x),r). 

Then u is a solution (always assumed BUG ) of 

d 

dtu + L (t, x, u, Du, D 2 u) = 2, (Du ■ Oi (x)) z\\ u (0, .) = u (.) 

if and only if u 1 (£, x) := u (£, ip (£, x)) is a solution of 

dt+L*>=0; u 1 (0, .) = «o (•) (8) 

Proof. Set y = ip t (x). When u is a classical sub-solution, it suffices to use the chain-rule and 
definition of F^ to see that 

w(t,x) = u (t, y) + Du (t, y) • V> t (x) = u (t, y) + Du (t, y) ■ V (y) z t 

< F (t, y, u (t, y) , Du (t, y) , D 2 u (t, y)) = F^ (t, x, w (t, x) , Dw (i, x) , D 2 w {t, x)) . 

The case when u is a viscosity sub-solution of (HJ) is not much harder: suppose that (i, x) is a 
maximum of w - £, where £, £ C 2 ([0,T] x K n ) and define tp e C 2 ([0,T] x K") by (p(t,y) = 
€ (*> ^t 1 (y)) ■ Set y = 4>t (x) so that 

F (i, y, u (i, y) , Dip (i, y) , D 2 y (t, y)) = F^ (i, x, w (£, x) , D£ (t, x) , D 2 £ (t, x)) . 

Obviously, (£, y) is a maximum of u — if, and since u is a viscosity sub-solution of (J6J we have 

<p (£, y) + Dip (£, y) V (y) z (t) < F (i, y, u (i, y) , Dip (t, y) , D 2 <p (i, y)) . 

On the other hand, £ (t, x) = p (t, ip t (x)) implies £ (t, x) = ip (£, y) + Dip (i, y) V (y) z (t) and putting 
things together we see that 

£ (l x) < F^ (t, x, w (i, x) , D£ (t, x) , D 2 £ (t, x)) 

which says precisely that w is a viscosity sub-solution of ([7]). Replacing maximum by minimum 
and < by > in the preceding argument, we see that if u is a super-solution of (|6]), then w is a 
super-solution of ([7]). 

Conversely, the same arguments show that if v is a viscosity sub- (resp. super-) solution for ([7|), 
then u(t,y) = w (t,ip~ (y)) is a sub- (resp. super-) solution for ([6]). ■ 

We prepare the next lemma by agreeing that for a sufficiently smooth function </> = </> (£, r, x) : 
[0,T]xlxl"4M wc shall write 

• d(j)(t 7 r 7 x) , dcf)(t,r,x) 



dt dr 

D^ : 1*^) and^ 2 --^ 20(t ' r ' a) 



dx% J l= l,...,n V dx*8x3 /iji=li ... in 



Lemma 5 (Outer transform) Let <j> = (p(t,r,x) G C 1 ' 2,2 and strictly increasing in r. (It follows 
that <j)(t,.,x) is an increasing diffeomorphism on the real line). Then u is a subsolution of dtu + 
F (t, x, r, Du, D 2 u) — 0, u (0, .) = uq (.) if and only if 

v (t, x) — <fT (t, u (t, x) , x) 

is a subsolution of dtV +^ F (t, x, r, Dv, D 2 vj = 0, v (0, .) = <\T (0, uq (x) , x) with 

«F(t,x,r,p,X) = ^ + ^F(t,x,tf>,D<f> + <f>'p, (9) 

<f>"p p + D<f> ®p + p® Defy' + D 2 (f> + <(>X) 

where <fi and all derivatives are evaluated at (t,r,x). The same statement holds if one replaces the 
word subsolution by supersolution throughout. 

Proof. (=>) Wc show the first implication, i.e. assume u is a subsolution of dfU + F = and set 

v (t, x) = (j)~ l (t, u (t, x) , x). By definition, (a,p, X) E V 2 ' + v (s, z) iff 

v (t, x) < v (s, z)+a (t — s)+p-(x — z)+ — (x — z) -X-(x — z)+o (\t — s\ + \x — z\ ) as (t, x) — > (s, z) . 
Since <f>(t,.,x) is increasing, 

and using a Taylor expansion on cf> in all three arguments we see that the rhs equals 

4> 0, v (s, z),z) + 8> „( SiZ)iZ (t - s) + (f)' sMSjZhz a (t- s) + 4>' s>v{SiZ)iZ p -{x-z) + ^' s ,v(s,z),z {x- zf -X -{x- z) 

+ D <Ps,v(s,z),z ■ i X - Z ) + 2 ( X ~ ^ ■ D V S ,„(«,*),* -(X-Z) 

+ {x-z) T -(D(<t ) ')) svMz ®p-{x-z) 
+ {x-z) T -p® {D(j))' sv{sz)z -{x-z) 

+ (x — z) ■ <j>" iV ( SiZ ) iZ p®p- (x — z) + o (\t — s\ + \x — z\ J as (s, z) — >• (t, x) 
Hence, 

\4>s,v(s,z),z + 4>s,v(s,z),z a i ^ >( t>s,v{s,z),z + ^S,t)(s,z),Z-P7 

4>",v(s,z),zP ®P + D {4>') StV{S:Z)tZ ®P+P® {D<l>)'s,v(8,z),z + D2( f>s,v(s,z),z + <P's,v(s,z),z X ) 

belongs to V 2;+ u (s, z) and since u is a subsolution this immediately shows 

4>s,v(s,z),z + ( t>(s,v(s,z),z) a + F (*> z i < t>(s,v(s,z),z)^'t>s,v(s,z),z + ^a,v{a,z),zPi 
4>",v(s,z),zP ®P + D (0') S ,„( S ,*),* ® P + P ® ( D 4>)'a,v(a,z),z + D *<i>8,v(<>,*),* + ^'s,v(s,z),z X ) < °- 



Dividing by <\> (> 0) shows that v is a subsolution of dtv + F^ = 0. 

(<==) Assume v is a subsolution of dfV +^ F = 0, ^F defined as in (|9]) for some F. Set u (t, x) :— 
4> (t, v (t, x) , x). By above argument we know that v is a subsolution of * (^F) (*> x j r jPi X) . For 
brevity write ip (t, ., x) = 0~ (i, ., ir) . Then 

^(^(t^r^.X) 

= -jf^ + —, ^ (t, X, 1pr t , r ,x)> D i>t,r,x + *P't,r,xP, 

Vt,r,x 4>t,r,x V 

<r,*P ® P + D (V»0 t r x ® P + P ® W)i,r,* + ^t.r,* + l^x* 



^t.r.a; 






Vt,r,x Vt.r,: 

T7 F (*. £> (*. V't.r,*' *) i £0*,^ r „x + ^Ut r *>* {-°^V,x + V't.r.xP} 

9t,i/) ttr>x ,x 

> ^^, r ,.,xP ® P + ^ (^)t,V- t , r ,„x ® P + P ® W)U.r...» + D2< f>t,i, ttria * 

+^ (ir ,.,x {^"r.zP ® P + ^ 0') i)riX <8> P + P <8> (£>V0U* + ^^W + <r,**}) 

Using several times equalities of the type (/ o f^ 1 ) = /*_i (/ _1 ) = «d cancels the terms involving 
4>, ip and their derivatives and we are left with F, i.e. 

t' 1 (*F) = F. 

This finishes the proof. ■ 

Corollary 6 (Transformation 2) Assume v = (r/x, . . . , Vd) C C°' 2 ([0,T] xl"). Assume 4> = 
4> (t, x, r) is determined by the ODE 

d 
<j) = (j) 2_j v i (*' x ) tf = 4> v (t, x) ■ it, (f> (0, x, r) = r. 

3=1 

Define L = L (t, x, r,p, X) by 

L= -Tr[A (t, x) ■ X]+b(t,x) -p + c (t, x, r) ; 

define also 

*L(t,x,r,X) = -Tr[A(t 1 x)-X]+' t 'b{t,x)-p+ 4 'c{t,x,r) where (10) 

2 
H(t,x)-p = b(t,x)-p T Tr[A(t,x)'D<t>'®p] 

*c(t,x,r) = -^Tr[A(t,x)-{D 2 cp)\+^b{t,x)-{Dcp) + ^ C {t,x,4>) 

where <j) and all its derivatives are evaluated at (t, r, x). Then 

dtw + L (t, x, r, Dw, D w) — w v (t, x) ■ z (t) = 
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if and only if v (t, x) = cf> (t, w (t, x) , x) satisfies 

d t v + L (t, x, r, Dv, D 2 v) = 0. 
Proof. Obviously, 



<(t,x,r) = rexp / VV(s,x)i;M 



This implies that <j>' — 4>/r and D<f> do not depend on r so that indeed ^b (t, x) defined above has 

no r dependence. Also note that cf>" — and <f>/4> — d- z = Ylj=x dj (t> x) %{• •"•* follows, for general 
F, that 

F(£,x,r,p,A) = rd-z + - 7 F(t,x,(j),D4> + (j)'p 1 

D(f> <g> p + p <g> L>0' + D 2 cj) + 0'A) 
and specializing to .F = L — wv ■ z, of the assumed (semi-) linear form, we see that 

*L = ---Tr [A (t,x)-(D<p' ®p + p(g> D<f>' + D 2 (f> + 0' X)] 

+ -jb (t, x) ■ (D4> + <p'p) + -jc (t, x, 4>) 

where </> and all derivatives are evaluated at (t, r, x) . Observe that ^L is again linear in X and p. 
It now suffices to collect the corresponding terms to obtain (fTU|) . ■ 
We shall need another (outer)transform to remove additive noise. 

Lemma 7 (Transformation 3) Set ip (t, x) = J Q g (s, x) ■ z s . Define 

L(t,X,r,p,X) = -Tr{A(t,x)-X]+b(t,x)-p + c(t 7 x,r); 
L a (t,x,r,p,X) = -Tr[A(t,x) ■ X]+b(t,x) -p + c a (t,x,r) 
with c a (t, x, r) = Tr [A (t, x) ■ D a (t, x)\ — b (t, x) ■ Da (t, x) + c(t,x,r — a (t, x)) 

Then v solves 

d t v + L (t, x, v, Dv, D 2 v) -g(t,x)-z (t) = 

if and only if v (t, x) = v {t, x) + (p {t, x) solves 

d t v + L a (t, x, v, Dv, D 2 v) = 0. 
Proof. Left to reader. ■ 

3.2 The full transformation 

As before, let 

L (t, x, r,p, X) := -Tr [A (t, x) X] + b(t,x) ■ p + c (t, x, r) 
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where A : [0, T] x R™ ->■ §", 6 : [0, T] x 1" -> R", / : [0, T] x R n x R -> R. Let us also define the 
following (linear, first order) differential operators, 

Mk(t,x,u,Du) = <Tk (%) ■ Du for k = 1, ... , di (11) 

Md x +k{t,X,u,Du) = uvk (t,x) for fc = 1, . . . , d 2 
M dl+d2+k (t,x,u,Du) = gk(t,x) for fc= 1, ...,d 3 . 

The combination of transformations 1,2 and 3 leads to the following 

Proposition 8 Lei z 1 G C 1 ([0,T],R dl ) , cr = (a u ...,a dl ) C C 6 2 (R n ,R") and denote the ODE 
flow ofdy = a (y) dz with ip, i.e. ip : [0, T] x R"^ R" satisfies 

i> (t, x) = o{ii (t, x)) z\, i\> (0, a;) = x G 1". (12) 

Further, let z 1 G C 1 ([0,T],R d2 ) and Zei i/ = {v x ,...,v dz ) be a collection o/C°' 2 ([0,T] xR",R) 
junctions and define </> = </>(£, r, x) as solution to the linear ODE 

4> = 4>v(t,iP t (x))%, (f>(0,r,x) =r el. (13) 

=d*(t,x) 

Further, let z 3 G C 1 ([0,T] ,R da ) and define ip (t,x) as the integral 

t 

4> n 4> 



where V (t >iC ) = -j^-^g (t, ^ (»)) 



?(*,&)= / V («,*)•#, (14) 

Jo 

1 
(/>'(*, a;) 

yl< toi, set z t := [z\,z 2 , zf) G R dl © R d2 © R d3 = R d . Then u is a viscosity solution of 

d t u + L (t,x,u,Du,D 2 uj = A (t,x, u,Du) z t , (15) 

u(0,x) = u (x), 

iff u (t, x) = (\T (t, u (t, ijj (t, x)) , x) + <p (t, x) is a viscosity solution of 

d t u + L(t,x,u,Du,D 2 u) = (16) 

u (0, x) = uo(x) 

where L = \j \L^y\ is obtained via transformations 1,2 and 3 (in the given order). 

Remark 9 Transformation 2 and 3 could have been performed in one step, by considering 

4> = <j>v* (t, x) ■ z 2 + g^ (t, x) -z\, <f> (t, r, x) | t=0 = r. 

Indeed, the usual variation of constants formula gives immediately 

(j){t,x)=rexp( [ ^(s,x)dz 2 \ + f eO^'^'V ( s , x) ■ dz 3 s 

and one easily sees that transformations 2 and 3 just separate the effects of this; with the benefit of 
keeping the algebra somewhat simpler (after all, we want explicit understandings of the transformed 
equations). 



3 Since is linear in r, there is no r dependence in its derivative <j>' . 
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Remark 10 Related to the last remark, generic noise of the form H (x, t, u) dz can be removed with 
this technique. The issue is that the transformed equations quickly falls beyond available theory, cf. 

m 

Proof. We first remove the terms driven by z 1 : to this end we apply transformation 1 with 
L (t, x, r,p, X) replaced by L — rv ■ z 2 — g ■ z 3 . The so-transformed solution, u 1 (t, x) = u (t, ipt ( x ))> 
satisfies the equation 

(d t + L*) u x -u x v (t, ^ (a;)) ■ z 2 - g (t, ^ (x)) ■ z 3 = 

=d*(t,a:) =c*(t,x) 

We then remove the terms driven by z 2 by applying transformation 2 with lr — c^ ■ z 3 . The so- 
transformed solution u 12 (t, x) = tfT x (t, u 1 (t, x) , x) satisfies the equation with operator (d t + & {ifi — g^ ■ i 3 )) 
i.e. 

d t v+ *(L^)v - 4/ -z 3 = 0. 



It now remains to apply transformation 3 to remove the remains terms driven by z 3 . The so- 
transformed solution is precisely u, as given in the statement of this proposition, and satisfies the 
equation 

(ft +[*(**)] J fi = o. 
The proof is now finished. ■ 

3.3 Rough transformation 

We need to understand transformations 1,2,3 when (z 1 , z 2 , z 3 ) becomes a rough path, say z. There 
is some "triadiagonal" structure: (fT2"j) can be solved as function of z 1 alone; 



dtp t (x) — a {ip t (x)) dz t with ip (x) = x. 



13"l) is tantamount to 



4> {t, r, x) = r exp 



v(s,ip s (x))dzl 



(17) 
(18) 



As for ip = <p (t, x), note that 



\/4> (t, r,x) — 4> (t, x) = exp 



v ( s > </> s 0*0) dz 2 s 



so that 



<p(t,x)= / (j>{s,x)g(s,ip s (x))dz 3 . 



(19) 



Lemma 11 Let z be a geometric p-rough path; that is, an element in C°' p ~ var ([0,T] , G^ (R d )). 
Let 7 > p > 1. Assume 

a = (a 1 ,...,o- dl )cLip^(W\R n ), 

v = (i/ 1 ,...,^ 2 )cIi^- 1 ([0,T]xR" ) R), 

9 = (gi,...,g d3 )cLtp^- 1 ([Q,T]xR n ,R). 
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Then ip,(j> and if depend (in local uniform sense) continuously on (z},z 2 ,z 3 \ in rough path sense. 
Under the stronger regularity assumption 7 > p+2; this also holds for the first and second derivatives 
(with respect to x) ofip,ip~ ,(j>,4> o,nd ip. In particular, we can define if),<j> and <p when (z 1 , z 2 , z 3 ) 
is replaced by a genuine geometric p-rough path z and write ?/> z , </> z , ip z to indicate this dependence. 

Proof. Given z one can build a "time-space" rough path (t, z) of (t, z 1 , z 2 , z 3 ) since the additionally 
needed iterated integrals against t are just Young integrals. Define 



W x 



The assumptions on a, v and g guarantee that W = {W lf ..., W A ) ■ R 1 + d ^+ 2d ^+ d 3 _>. L (R 1+d , Ri+<*i+2<*2+<*3) 
is Lip 1 . Hence, the "full" RDE (parametrized by x G R n and r e M) 

/ 1 \ 

cr(V>) 

d(t,z)= r.0.i/(*,^) d(t,z) 

-<j>.v(t,ip) 

V hM) I 

has a unique solution (with obvious initial condition that the flows ip,4>,ip evaluated at t — are 
the identity maps). Further, every additional degree of Lipschitz regularity allows for one further 
degree of differentiability of the solution flow with corresponding stability in rough path sense. ■ 

Remark 12 The regularity assumptions with respect to t could be relaxed; see the forthcoming work 
by Caruana and Gurko. 

4 Semirelaxed limits and rough PDEs 

The goal is to understand 





( l \ 




( l \ 




V> 




v> 


d 


<t> 


= w 


<P 














U/ 




\<p J 



dl 



d t u + L (t, x, u, Du, D u) = y [pi (x) ■ Du) z t ' 1 + u\ Vj (t, x) z t ' J + > gt (t. 



X) z 



?,:k 



i=l 



fc=l 



in the case when (z 1 , z 2 , z 3 ) becomes a rough path. To this end we first need quantitative assump- 
tions on L; we shall assume that 



L(t,x,r,p,X) 



-Tr 



a (t, x) a (t, x) X +b(t,x) ■ p + f (t,x,r) 



(20) 



where a (t, x) : [0, T] xlMK" x " and b (t, x) : [0, T]xKM M™ are bounded, continuous in t and 
Lipschitz continuous in x, uniformly in £ e [0, T]. We also assume that / : [0,T] x M™ x M — > R is 
continuous, bounded whenever r remains bounded, and with a lower Lipschitz bound, i.e. 3C < 
s.t. 

/ (t, x, r) - f (t, x,s)>C(r- s) for all r > s, (t, x) e [0, T] x W 1 . (21) 
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This will guarantee that comparison holds for dt + L; see the appendix for some details. Let us now 
replace the (smooth) driving signals of the earlier sections by a {d\ + di + d 3 )-dimensional driving 
signal z £ and impose convergence to a genuine geometric p-rough path z. That is, in the notation 
of Chapter 14] 

z g C* ^ var ([0,T] ,G [p] (E. dl+d2+d3 )\ . 

We shall write r ^' L ,4) z ,a z for the objects objects (solution of rough differential equations and in- 
tegrals) built upon z, as discussed in the last section (lemma fTTj) . We also write ijj e ,(f) E ,a e when 
the driving signal is z e . Recall from ([TT|) that A = (Ai, . . . , Ad) is a collection of linear, first order 
differential operators. We have 



Theorem 13 Let p > 1. Let L be of the above form 120\) with assumptions on <j, 6, / as details 
above. Assume furthermore that the coefficients of A — (Ai, . . . , A e i 1 +d 2 +d 3 ) have Lip 1 -regularity 
where 7 > p + 2. Let uq £ BUC and z a geometric p-rough path. Then there exists a unique 
u = u z e BUC ([0, T] x R n ), write formall$3, 

du + L (t, x, u, Du, D 2 u\ dt = A (t, x, u, Du) dz, u (0, •) = uq 

such that for any smooth sequence z e — > z in p-rough path sense, u e , viscosity solutions (always 
assumed BUC) of the classical PDE problems 

if + L (t, x, u £ , Du 6 , D 2 u e ) = A (t, x, u e , Du e ) z e , u (0, •) = u , 

converge locally uniformly to u z . Moreover, we have the contraction property 






:x :J 



(C given bv \21)) and continuity of the solution map (z,« ) •->• u z from 

C p-var A Qi T ] ^ G [p] £ R d)\ x BUC ( R nj _^ BUC Qq^ j,] x R n) _ 

Proof. We use the same technique of "rough semi-relaxed limits" as in [CFO10]: the key remark 
being that 

'<r(L*'X\ ^_> [^ (l^)} r 

holds locally uniformly, as function of (t,x,r,p, X). Secondly, the classical solutions u e can be 
transformed to 

V e (t, x) = {^Y 1 (t, u £ (t, iP £ (t, x)) ,x)+a e (t, x) ■ 

it is a simple matter, with comparison for instance to see that, for suitable C, 

sup \v e {t,x)\ < (l + |u | oo )e c ' T ; 
ee(o,i] 
te[o,T] 

x£S. n 

this in turn implies (thanks to the uniform control on a e ,<jf,^i e as e — > 0, using the rough path 
representations discussed in section |3~3|) that u e (t,x) remains locally uniform bounded (as e — > 0). 



10 The intrinsic meaning of this "rough" PDE is discussed in definition 1141 below. 
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The reader may wonder if u is the solution in a sense beyond the " formal" equation 

du + L (t, x, u, Du, D u) dt = A (t, x, u, Du) dz, u (0, •) = uq. 

Inspired by the definition given by Lions-Souganidis in [9] we give 

Definition 14 u is a solution to the above rough partial differential equation if and only if 

dv + L (t, x, v, Dv, D 2 v) = 0, v (0, •) = Uo 

in viscosity sense where 

L = \ r (l^) . 

V / J a z 

Corollary 15 There exists a unique solution to the rough PDE. 

Proof. Existence is clear from theorem [T3] Uniqueness is inherited from uniqueness to the Cauchy 
problem for (dt+L\ =0 which follows from comparison. ■ 

5 Appendix: comparison for parabolic equations 

Let u £ BUC ([0, T] x R") be a subsolution to d t + F; that is, d t u + F (t, x, u, Du, D 2 u) < if u is 
smooth and with the usual viscosity definition otherwise. Similarly, let v G BUC ([0,T] x R") be a 
supersolution. 

Theorem 16 Assume condition (3.14) of the User's Guide |2J/ ; uniformly in t, together with uni- 
form continuity of F — F (t,x,r,p, X) whenever r,p, X remain bounded. Assume also a (weak form 
of) properness: there exists C such that 

F(t,x,r,p,X)-F(t,x,s,p,X)>C(r-s) Vr < s, (22) 

and for all t G [0, T] and all x,p,X. Then comparison holds. That is, 

u (0, •) - v (0, •) => u < v on [0, T] x R". 

Proof. It is easy to see that u = e~ u is a subsolution to a problem which is proper in the usual 
sense; that is (|2"2"j) holds with C = which is tantamount to require that F is non-decreasing in r. 
The standard arguments (e.g. [3] or appendix of 2]) then apply with minimal adaptations. ■ 

Corollary 17 Under the assumptions of the theorem above let u, v be two solutions, with initial 
data uq,vq respectively. Then 

\ u - w loo;R"x[0,T] < e ° T Wo - WOloojH- 

with C being the constant from l\22t . 

Proof. Use again the transformation u = e~ ct u,v — e~ ct v. Then v + \uq — woIoo-r™ is a superso- 
lution of a problem to which standard comparison arguments apply; hence, 



U < V + \U - VoloojRn • 

Applying the same reasoning to u+ \uq — ^oloo-jj" an d finally undoing the transformation gives the 
result. ■ 
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